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XEFT is a low-energy effective field theory for charm mesons and pions that provides a systemat¬ 
ically improvable description of the X (3872) resonance. A Galilean-invariant formulation of XEFT 
is introduced to exploit the fact that mass is very nearly conserved in the transition —>■ 

The transitions —>■ D°tt° and X —^ are described explicitly in XEFT. The effects of 

the decay D*° —>■ D °7 and of short-distance decay modes of the X(3872), such as , can 

be taken into account by using complex on-shell renormalization schemes for the D*° propagator 
and for the D*°D° propagator in which the positions of their complex poles are specified. Galilean- 
invariant XEFT is used to calculate the D*^D° scattering length to next-to-leading order. Galilean 
invariance ensures the cancellation of ultraviolet divergences without the need for truncating an 
expansion in powers of the ratio of the pion and charm meson masses. 

PACS numbers: 14.40.Rt, 14.40.Lb 


I. INTRODUCTION 

The surprising discovery of the X(3872) by the Belle Collaboration in 2003 [1] marked the beginning of a renaissance 
in quarkoninm spectroscopy [2]. Dozens of new mesons whose constituents include a heavy quark and antiquark and 
with mass above the open-heavy flavor threshold have been observed. They are collectively referred to as XYZ mesons. 
Some of the XYZ mesons are electrically charged and therefore must be tetraquark mesons whose constituents also 
include a light quark and antiquark. The pattern of the observed XYZ mesons remains unexplained. In particular, 
the relation between the X(3872) and the other XYZ mesons is still not understood. 

The discovery decay mode of the Ar(3872) implies that its constituents must include a charm quark and 

antiquark. However the Al(3872) has properties inconsistent with conventional charmonium, including comparable 
branching fractions into decay modes with isospin 0 and isospin 1. The quantum numbers of the X(3872) were 
finally established by the LHCb Collaboration in 2013 to be l+“'' [3]. Its mass is extremely close to the threshold 
for the pair of charm mesons By combining precise measurements of the mass Mx of the Ar(3872) in the 

J/V'TT+TT” channel with precise measurements of the masses M* and M of the D*^ and D^, the difference 6x between 
the threshold and the mass has been determined to be 

Sx = iM„ + M)-Mx = 0.11 ± 0.23 MeV. (1) 

The quantum numbers I’*"'*" of the X(3872) imply that it has an S-wave coupling to the charm meson pairs D*°D° 
and Given the small value of Sx, the universality of near-threshold S-wave resonances implies that the X(3872) 

must be a bound state (if 6x > 0) or a virtual state (if Sx < 0) whose constituents are the C = + superposition 
D*of)0 _|_ observation of X(3872) in hadron collisions strongly suggests that it is a bound state (like 

the deuteron) rather than a virtual state (like the dineutron). One universal property of S-wave near-threshold bound 
states is that the mean separation {r)x of the constituents is determined by the binding energy Sx- {r)x = 
where is the reduced mass of [5]. Given the binding energy Sx = O.ll+o'^a MeV, the mean separation of the 

charm mesons is predicted to be 6.8^“g fm. Thus the size of the X(3872) is comparable to that of the largest nuclei. 

The universal properties of S-wave near-threshold bound states imply that most of the probability of the X(3872) 
is in a molecular component consisting of well-separated charm mesons D*^D^. (From now on, the equally probable 
D^D*^ component will usually not be mentioned explicitly.) At short distances, the wavefunction of the A(3872) 
can have other components with smaller probabilities. One possible component is the 1“'"+ P-wave charmonium state 
Xci(2T’), whose constituents are cc. Another possibility is an isospin-0 tetraquark with constituents cqcq, where g is a 
light u 01 d quark. It could be a compact tetraquark, it could have substructure consisting of the diquark clusters cq 
and cq, or it could have substructure consisting of color-singlet clusters, such as the pair of mesons J/if uj. There are 
also hexaquark components of the wavefunction with constituents cqqcqq. One such component that is particularly 
important is , because the constituent D*^ in the dominant component of the wavefunction can decay into 

. In considering the various possible components of the A(3872) wavefunction, it is essential to take into account 
their couplings to D*^D^ and the resonant interactions between the charm mesons. 
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The interplay between the D*^D^, , and other components of the wavefunction can be treated systematically 

using an effective field theory in which the and 7r° are explicit degrees of freedom. Such an effective field 

theory has been developed by Fleming, Kusunoki, Mehen, and van Kolck and named XEFT [6]. In its simplest form, 
XEFT is a nonrelativistic field theory for and 1?°, their antiparticles, and 7r°. It is straightforward to extend 
XEFT to include the charged charm mesons and the charged pions. In XEFT, the interactions between and 
must be treated nonperturbatively in order to generate a bound state that can be identified with the X(3872). 
Fleming et al. showed that pion-exchange interactions can be treated perturbatively along with the range corrections 
to D*D interactions [6]. We denote the masses of and 7r° by M*, M, and m, respectively. We denote the 

difference between the mass and the sum of the and 7r° masses by 5: 

(5 = M* - (M + to) « 7.14±0.07 MeV. (2) 

The power counting of XEFT is defined by taking external momenta, the binding momentum scale ^/M6x of the 
X(3872), and the momentum scale V mS of the pion to all be low momentum scales of order Q. The low energy 
scales include the binding energy dx of the X(3872), the kinetic energy scale S for a pion, and the kinetic energy 
scale mS/M for a charm meson. The high momentum or energy scales include m, M, and dTr/^r, where is the pion 
decay constant. Amplitudes are calculated as systematic expansions in powers of Q divided by a high momentum or 
energy scale. In the original paper on XEFT by Fleming et ah, the momentum distributions in the decay of A(3872) 
to were calculated to next-to-leading order (NLO) in the XEFT power counting and to leading order in an 

expansion in powers of m/M [6]. 

There have been a number of subsequent applications of XEFT. Fleming and Mehen applied XEFT at leading 
order (LO) to decays of the X = A(3872) into the P-wave charmonium state XcJ plus one or two pions [7, 8]. Mehen 
and Springer applied XEFT at LO to the radiative decays X —>• tp{2S)x and ^’(4040) —>■ Ay [9]. Margaryan and 
Springer applied XEFT at LO to the decay ^/>(4160) —>■ Ay [10]. Recently, Jansen, Hammer and Jia used XEFT to 
determine the dependence of the binding energy of A on the light quark masses [11]. They also calculated the D*^D^ 
scattering length to NLO and to leading order in m/M. The applications of XEFT are not limited to the DDtt sector, 
which also includes D*D and DD* states. Canham, Hammer, and Springer pointed out that XEFT can be applied 
to the DDDtt and DDDtttt sectors [12]. They calculated the S-wave phase shifts for low-energy scattering of D or 
D* from A at LO. Braaten, Hammer, and Mehen pointed out that XEFT can be applied to the DDtttt sector, which 
also includes D*D*, D*Dt:, and DD*tt states [13]. They calculated the low-energy cross sections for elastic tt+A 
scattering and for the break-up reaction tt+A —>■ at LO. 

We will refer to the formulation of XEFT presented in Ref. [6] as original XEFT. There are various problems 
with original XEFT that present obstacles to accurate quantitative predictions. Some of these problems are due 
to its formulation as a nonrelativistic field theory of the charm mesons and 7r° that is not Galilean invariant. One 
consequence of the lack of Galilean invariance is that a preferred frame (such as the center-of-momentum frame) must 
be specified either explicitly or implicitly in any calculation. Another consequence of the lack of Galilean invariance 
is that ultraviolet divergences are much less constrained. The explicit counterterms in the Lagrangian for original 
XEFT in Ref. [6] are sufficient to eliminate ultraviolet divergences in NLO calculations only if results are expanded 
in powers of m/M and then truncated at a sufficiently low order. This truncation provides a limit on the accuracy. 
Although the mass ratio m/M = 0.072 is very small, the expansion is actually in powers of the square root of the mass 
ratio, which is 0.27. The alternative of adding additional counterterms to the Lagrangian to cancel the ultraviolet 
divergences would introduce additional parameters that would have to be determined phenomenologically. Another 
problem with original XEFT is that the renormalization scheme made it difficult to take into account some decays into 
final states with momenta too large to be treated explicitly in XEFT. These decays are D*^ —>■ iJ°y, which accounts 
for about a third of the full width of the D*^, and all decay modes of the A(3872) other than . 

Alhakami and Birse have recently proposed an alternative power counting for XEFT [14]. In their power counting, 
m and M* — M are treated as small energy scales of order Q, and 5 is treated as a tiny energy scale of order Q'^/M. 
The ratios 6/m = 0.051 and m/M = 0.072 are both of order Q/M. This power counting scheme was not implemented 
at the lagrangian level, but it provides an organizing principle for simplifying Feynman diagrams. The power counting 
of Ref. [14] makes the expansion in powers of m/M systematic. The formulation in Ref. [14] is convenient for matching 
onto heavy hadron chiral perturbation theory, which has been used in several applications of XEFT [7-10]. It does 
not address the problems of frame dependence or of additional ultraviolet divergences. With this new power counting, 
the decay into is within the domain of applicability of XEFT. However there are important decay modes of 
the A(3872), such as J/ 1 /, that remain outside the domain of applicability of XEFT. 

In this paper, we present a new formulation of XEFT that removes many of the obstacles to accurate quantitative 
calculations. The new formulation is an effective field theory with a Galilean symmetry that is motivated by the fact 
that mass is very nearly conserved in the transition —>■ The Galilean symmetry solves the problem of 

frame dependence and it dramatically simplifies ultraviolet divergences. To take into account the decay D*^ —>■ 
and decay modes of the A(3872) other than a new renormalization scheme is introduced that is expressed 
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in terms of the complex energies of D*^ and X(3872). The new Galilean-invariant formulation of XEFT is illustrated 
by a calculation of the D*^D^ scattering length to NLO in the XEFT power counting. The cancellation of ultraviolet 
divergences for arbitrary values of m/M is verified explicitly. 


II. GALILEAN-INVARIANT XEFT 

In this section, we introduce a new formulation of XEFT as a Galilean-invariant field theory. We describe the 
changes in the Lagrangian for the original XEFT defined in Ref. [6] that are required for Galilean invariance. We 
explain how the effects of the decay D*^ —>■ can be taken into account through the complex rest energy of the 

D*^. We then write down the next-to-leading order Lagrangian for Galilean-invariant XEFT, and give its Feynman 
rules. 


A. Galilean Invariance 


A unique feature of the decay D*^ —>■ is that mass is very nearly conserved: the sum of the masses of the 

and 7r° is only about 3.5% lower than the mass of the Galilean invariance is a possible space-time symmetry of 
a nonrelativistic theory that requires exact mass conservation [15]. The mass that must be conserved is the kinetic 
mass, which is the mass that appears in the denominator of the kinetic energy. The approximate conservation of 
mass in the transition D*^ —>■ strongly motivates a Galilean-invariant formulation of XEFT. The Lagrangian for 

original XEFT, including all terms required to calculate the decay of A'(3872) into to next-to-leading order 

(NLO), was written down in Ref. [6]. We will describe how the terms in this Lagrangian must be modified to make 
them Galilean invariant. The NLO Lagrangian for Galilean-invariant XEFT will be written down in Section IIC. 

We choose the kinetic masses of the and 7r° to be their physical masses M and m, respectively. Conservation of 
kinetic mass then requires the kinetic mass of D*^ to be M -|- m. The difference S between the mass of the D*^ and 
its kinetic mass must be taken into account through its rest energy. Kinetic mass conservation requires the following 
change in the Lagrangian for original XEFT defined in Ref. [6]: 

• In the kinetic term ■ V‘^D/{2mD*) for the D*'^, its mass mo* must be replaced by M -|- m. 


A similar change must be made in the kinetic term for the It will be convenient to introduce the reduced kinetic 
mass /i for and the reduced mass for 




Mtt 


M {M + m) 


2M + m 

mM 


= 965.0 MeV, 


M + m= 125.87 MeV 


(3a) 

(3b) 


The ratio of these reduced masses is 


r= — = 0.1304. 


( 4 ) 


Galilean invariance requires interaction terms to be invariant under Galilean boosts, in which the momenta of tt*^, 
ZZ°, and ZZ*° are boosted by a common velocity vector v multiplied by their kinetic masses m, M, and M -\- m, 
respectively. Galilean invariance requires three changes in the interaction terms in the Lagrangian for original XEFT 
that given in Ref. [6]: 


• In the pion interaction term ■ D'Vtt, the operator V between D and tt 
m ^)/{M + m). 


should be replaced by (M^ — 


• In the D*^D^ interaction term {DD^ ■ D{^YD, the operator (^Y between D and D should be replaced 
by 4(M^ — (M -I- m)^)^/(2M -|- niY- 




In the interaction term {DDY ■ ZZZZVtt that describes the transition of Z?°Z?°7r° 
between DD and tt should be replaced by (2M^ — m ^)/(2M + m). 


to the operator V 
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Similar changes must be made in the hermitian conjugates and charge conjugates of these interaction terms. The 
three modified interaction terms described above reduce to those of original XEFT in the limit m/M —>■ 0. The 
modified interaction terms ensure the invariance of amplitudes under Galilean boosts. 

Galilean invariance also strongly constraints ultraviolet divergences. The changes above specify all terms in the 
Lagrangian for Galilean-invariant XEFT that are required to calculate the decay of X(3872) into to NLO in 

the XEFT power counting. The results should be independent of the ultraviolet cutoff to all orders in m/M. Without 
Galilean invariance, there are three independent D*^D^ interaction terms and two independent —>■ D*^D^ 

interaction terms. The Lagrangian for original XEFT defined in Ref. [6] includes only one interaction term of each 
kind. Cutoff independent results for the decay rate of X(3872) into at NLO were obtained by also truncating 

the expansion in powers of m/M at leading order. 

Galilean-invariant XEFT can be extended to include charged charm mesons and charged pions. The tt’’' must have 
the same kinetic mass m as 7r°. The must have the same kinetic mass M as 77°. The must have the same 
kinetic mass M -|- m as I?*°. The difference between the mass and kinetic mass of a particle must be taken into 
account through its rest energy. 

In a nonrelativistic effective field theory for charm mesons, one can impose a phase symmetry that guarantees 
the separate conservation of the number Nc of charm quarks and the number Ng of charm antiquarks. These quark 
numbers can be expressed in terms of meson numbers: 

Nc = Njyto -\-N£)0, ( 5 &) 

Nc = Nq,o+Nqo. (5b) 

If one also considers charged charm mesons, the charm quark number Nc is the sum of the numbers of D*^, D*^, 
11°, and . In Galilean-invariant XEFT, the exact conservation of kinetic mass in the transitions 77*° O 77°7r° and 
77*° o 77°7r° provides motivation for introducing an additional phase symmetry that guarantees the conservation of 
the pion number defined by 


Nt^ — N.^q -\- N -j- N. (6) 

The name pion number is appropriate since the D* can be interpreted as a P-wave Dtt resonance. In Galilean-invariant 
XEFT with charged charm mesons and pions, the pion number is the sum of the numbers of 7r°, 7r+, 7r“, 77*°, 77*°, 
77*+, and 77*". 


B. Choice of Rest Energies 

In a nonrelativistic field theory, deviations from mass conservation in a reaction are taken into account through 
the rest energies of the particles involved. The rest energy of a particle can be chosen to be its physical mass. In a 
Galilean invariant theory, since the kinetic mass is conserved, the rest energy of a particle can equally well chosen as 
the difference between its physical mass and its kinetic mass. If there are linear combinations of the particle numbers 
that are conserved, field redefinitions can be used to set the rest energies of some of the particles to 0. In XEFT, 
conservation of the charm quark number Nc allows the rest energies of 77° and 77*° to be set to 0 and M* — M, 
respectively. Conservation of both Nc and the pion number N^^ allows the rest energies of 7r°, 77°, and 77*° to be set 
to 0, 0, and 6 = M* — M — m, respectively. 

The rest energy of a particle can also be used to take into account its partial width into decay modes that cannot 
be described explicitly in the effective field theory. Such a partial width is taken into account through a negative 
imaginary term in the rest energy. The width of 77*° is more than 5 orders of magnitude larger than the widths of 
77° and 7r°, so the widths of 77° and 7r° can be completely ignored. The width of 77*° come from its decay into 77°7r°, 
which can be described explicitly in XEFT, and from its decay into 77°7, whose momenta are too large to be described 
explicitly. The effect of the 77°7 decay mode can be taken into account by including a term —7r*o,7/2 in the 77*° 
rest energy, where r*o ,7 is the partial width of 77*° into 77°7, whose branching fraction is approximately 38%. If we 
use the complex on-shell renormalization scheme for the 77*° propagator in which the position of the complex pole is 
specified, the 77*° rest energy also includes the term —7r*o,7i-/2, where L^o.tt is the partial width of 77*° into 77°7r°, 
whose branching fraction is approximately 62%. In this scheme, the complete imaginary part of the 77*° rest energy 
is —7r*o/2, where r,o = E^o.tt + r*o ,7 is the full width of the 77*°. 

In original XEFT, the rest energies of 7r°, 77°, and 77*° were chosen to be —S, 0, and 0, respectively. Setting the 
rest energy of 77*° to 0 is inconvenient if we wish to take into account the partial width of 77*° into 77°7. It is more 
convenient to choose the rest energies of 7r° and 77° to be 0. The rest energy of 77*° is then S — zr,o/2. This requires 
the following change in the Lagrangian for original XEFT: 
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• The rest-energy term (57r^7r for the 7r° should be replaced by the rest-energy term —(5 — zr,o/2)-D^ • D for the 

D*o 

There is a similar rest-energy term for the D*^. 

If Galilean-invariant XEFT is extended to include charged charm mesons and charged pions, mass differences and 
nonnegligible decay modes outside the domain of validity of the effective field theory can be taken into account through 
the rest energies of the particles. The 7r+ — 7r° mass difference can be taken into account through the rest energy of 
7r+. The D'^ — mass difference can be taken into account through the rest energy of . The difference between 
the D*^ — mass difference can be taken into account through the rest energy of The tiny partial width 

for the decay of D*^ into D+y can also be taken into account through the rest energy of 

C. NLO Lagrangian 

We can now write down the complete NLO Lagrangian for Galilean-invariant XEFT for the neutral charm mesons 
and the 7r°. It includes only those terms required to calculate the elastic scattering amplitude for to NLO in 

the XEFT power counting. The field for the 7r° is denoted by tt. The fields for the and are denoted by D and 
D. The fields for the D*^ and D*^ are denoted by D and D. The Lagrangian is the sum of kinetic terms, interaction 
terms, and counterterms. 

There are kinetic terms for the 7r°, I?°, , and I?*°, and The kinetic terms for the 7r°, I?°, and D*^ are 


£^o = (z9o + V^/(2m)) TT, (7a) 

Ldo = D^idQ + V^/{2M))D, (7b) 

Cd-o = • (z9o + vV(2(M-hTO))) D-((I-zr*o/2)r>'^ • r>. (7c) 

The kinetic terms for the and D*^ can be obtained from those for and D*^ by replacing D hy D and D by D. 

The interaction terms consist of pion interactions, contact interactions, and interactions. The pion interaction 
terms for the transitions D*^ O are 

JO-D‘O^DO.0 = f \d^ ■ iD[M^ - mV]7r) + (7?[M^ - mV]7r)t • d] . (8) 

lyTTiji^ylVl +m) I 


The pion interaction terms for the transitions e-?- D^tt^ are obtained by replacing D and D hy D and D. The 
contact and interactions are the same in the —>■ D*^D^, D^D*^ o D*^D^, and D^D*^ —>■ D^D*^ channels. 

The interaction terms in the —>■ D*^D^ channel are 

£d.odo^d-d" = ■ (6r>) 

T{D[M^-{M + -{DD)]. (9) 

The interaction terms in the other three channels are obtained by replacing D and D hy D and D in the appropriate 
places. The coupling constants g in Eq. (8) and Cq and C 2 in Eq. (9) are essentially the same as those in the 
Lagrangian for original XEFT defined in Ref. [6]. 

To cancel ultraviolet divergences in Green functions at NLO, it is also necessary to include counterterms in the 
Lagrangian. The counterterms are the same in the D*^D^ —>■ D*^D^, D^D*^ O D*^D^, and D^D*^ —>■ D^D*^ 
channels. The counterterms in the D*^D^ —>■ D*'^D^ channel are 

-^counterterm ^ (DD)^ ■{DD)-^iDDy-idtiDD). (10) 

The counterterm with coefficient SDq is not needed to calculate on-shell quantities, because it can be reduced to the 
other counterterm by using field redefinitions. It was omitted in Ref. [6], because the momentum distribution for 
in the decay of the X(3872) is an on-shell quantity. 
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FIG. 1: The propagators for tt®, D^, and are represented by dashed, solid, and solid+dashed lines, respectively. The 
Feynman rules for these propagators are given in Eqs. (11). The propagators for and D*^ look like those for and D*^ 
with the arrows reversed. 



FIG. 2: The contact interaction vertex for D*°D^ —>■ D*^D^. The Feynman rule for this vertex is given in Eq. (13). 


D. Feynman Rules 


The terms in the NLO Lagrangian for Galilean-invariant XEFT with neutral charm mesons and the tt^ are given in 
Eqs. (7), (8), (9), and (10). We will give the Feynman rules for this Lagrangian. The Feynman rules that are actually 
used for calculations beyond leading order will be enclosed in boxes. 

In XEFT, the charm quark and antiquark numbers W and iVg are conserved. In Galilean-invariant XEFT, the pion 
number iV^ is also conserved. These conservation laws can all be made manifest in the Feynman rules by appropriate 
notation for the propagators. We use a dashed line for the pion propagator, a solid line for the D and D propagators, 
and a double line consisting of a solid and a dashed line for the D* and D* propagators. In the propagators for the 
mesons D* and D that contain a charm quark, the solid line has a forward arrow. In the propagators for the mesons 
D* and D that contain a charm antiquark, the solid line has a backward arrow. It is sometimes convenient to omit 
the arrows on internal lines of diagrams and to use the convention that there is an implied sum over the possible 
directions of the omitted arrows. 

The propagators for the 7r°, 11°, and D*^ are illustrated in Fig. I. The Feynman rules for the propagators of the 
7r°, the 11° or 1?°, and the D*^ or D*^ are 


(lla) 

(llb) 

(llc) 

where po and p are the energy and momentum of the particle. The propagator of the D*^ is diagonal in its vector 
indices i and j, and its complex rest energy is 



E^=S-ir^o/2. ( 12 ) 

Since F* has a negative imaginary part, an explicit ie prescription is unnecessary in the D*^ propagator. 

At LO in the XEFT power counting, the only interaction that is required is a contact interaction for in 

the C = + channel. The Feynman rule for the vertex is the same for the —>■ and 

£)Ojj*o DOf)*o contact interactions: 


{-iCo/2)S^E 


( 13 ) 
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FIG. 3: The vertices for the pionic transitions D°tt° and ^ D*^. The Feynman rule for the D*° —>■ D°7r° vertex 

is given in Eq. (14). 




FIG. 4: The vertex and the counterterm vertex for D*^D° —>■ D*°D°. The Feynman rules for these vertices are given in 
Eqs. (15) and (16), espectively. 


This vertex for D*^D^ —>■ D*^D^ is illustrated in Fig. 2. The Cq interaction must be treated nonperturbatively in 
XEFT. The set of subdiagrams consisting of an arbitrary number of successive Cq interactions can be summed to all 
orders analytically. The Feynman rule for the resulting effective interaction is given later in Eq. (44). In calculations 
beyond leading order in XEFT, that effective vertex replaces the one in Eq. (13). That is why the Feynman rule in 
Eq. (13) is not enclosed in a box. 

One of the interactions beyond LO in the XEFT power counting is the D* o Dtt transition. The Feynman rules 
for the D*^ —>■ D^tt^ and —>■ D^tt^ vertices in Galilean-invariant XEFT are 
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{Mq — mpY 


M + m ’ 


(14) 


where q and p are the momenta of the outgoing 7r° and charm meson, respectively. The vertices for D*^ —>■ and 

—>■ D*^ are illustrated in Fig. 3. In the center-of-momentum frame defined by p -f q = 0, the momentum- 

dependent factor in Eq. (14) reduces to q*, which is the momentum-dependent factor in all frames in original XEFT. 
The Feynman rules for the —>■ and Z)°7r° —>■ vertices differ from Eq. (14) by an overall minus sign if 

q and p are the momenta of the incoming 7r° and charm meson, respectively. A convenient way to implement the 
XEFT power counting is to assign orders in g to the coupling constants of all other interaction terms. A complete 
calculation then requires calculating all diagrams to a given order in g. 

The other interaction at NLO in the XEFT power counting is the interaction for in the C = + channel. 

The Feynman rule for the vertex is the same for the D*^D^ —>■ D*^D^, D^D*^ O D*^D^, and D^D*^ —>■ D^D*^ 






interactions. In Galilean-invariant XEFT, the Feynman rule is 


/ ((^ + ^)P - ^P*f + ((^ + Hp' - ^ - 

-(2M + m)2-^ ’ 


(15) 


where p and p* are the momenta of the incoming spin-0 and spin-1 charm mesons and p' and p(. are the momenta 
of the outgoing mesons. The vertex for D*^D^ —>■ D*^D^ is illustrated in Fig. 4. In the D*^D^ center-of-momentum 
frame dehned by p* -f p = p(, -I- p' = 0, the momentum-dependent factor in Fq. (15) reduces to p^ -I- p'^. In 
this frame, the Feynman rule coincides with that in original XFFT, for which the momentum-dependent factor is 
[(p — p*)^ -I- (p' — p(,)^]/4. The coupling constant C 2 in Fq. (15) can be assigned order g'^. 

In addition to these interaction vertices, it is also necessary to include a counterterm vertex for in the C = + 

channel. The vertex is the same for the D*^D^ —>■ D*^D^, D^D*^ O D*^D^, and D^D*^ —>■ D^D*^ interactions: 


{-i[5Co + 5D^E]/2)5^\ 


(16) 


where E is the total energy of the pair of charm mesons. The counterterm vertex for D*^D^ —>■ D*^D^ is illustrated 
in Fig. 4. The SDq term is not needed to calculate on-shell quantities at NLO, but it is needed to cancel ultraviolet 
divergences in off-shell Green functions. The counterterm coefficients SCq and SDq in Fq. (16) can be assigned order 

In the original paper on XFFT, an additional interaction term that produces a transition of D*^D^ to was 

written down explicitly [6]. The Feynman rule for its vertex in Galilean-invariant XFFT is 


Bi {2Mq-m{p + p)y 
2y^ 2M + m ’ 

where q, p, and p are the outgoing momenta of the pion and the two spin-0 charm mesons. In the center-of- 

momentum frame defined by p -I- p -f q = 0, the momentum-dependent factor in Fq. (17) reduces to g*, which is the 
momentum-dependent factor in all frames in original XFFT. The coupling constant Bi in Eq. (17) can be assigned 
order . This interaction term was needed in Ref. [6] to calculate the decay of X(3872) into to NLO in the 

XEFT power counting or, equivalently, to relative order g^. 

In the D*^ propagator in Eq. (11c), a wavefunction renormalization factor Z can be inserted in the numerator. 
That factor can be made completely arbitrary by inserting canceling factors of into the interaction vertices for 

each D*'^ or Z)*° line. The number of such factors in the interaction vertices in Eqs. (13), (14), (15), and (16) would 
be 2, 1, 2, and 2, respectively. These factors can be absorbed into the coupling constants. 

If we use dimensional regularization with d spatial dimensions, it is useful to introduce a renormalization scale A 
to keep the dimensions of coupling constants the same as in the physical dimension d = 3. The interaction vertices 
in Eqs. (13), (14), (15), and (16) would be multiplied by A raised to the powers 3 — d, (3 — d)/2, 3 — d, and 3 — d, 
respectively. In a Green function, the net effect of these powers of A is a factor of A^""^ for every loop integral and 
an overall multiplicative factor of A^”'^ raised to some power. The factor of A^”"^ associated with a loop integral 
can be absorbed into its integration measure. If the Green function is made finite by renormalization, the overall 
multiplicative factor of A^”'^ raised to some power can be simply be ignored, because it is equal to 1 in the physical 
dimension d = 3. 


III. RENORMALIZATION PRESCRIPTIONS 

In this section, we introduce renormalization prescriptions that can be used to take into account the effects of 
transitions to states in which the momenta are too large to be described explicitly in XEFT. We also calculate the 
transition amplitude for D*^D^ —>■ D*^D^ at leading order in XEFT. 


A. D* 


Dn Coupling Constant 


One of the coupling constants in XEFT beyond leading order is the D* — Dn coupling constant g/\/2fT^. It could 
in principle be determined from the partial width of into 


= r[D*° ^ D° 


”1 ^ (34) 


(18) 
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FIG. 5: The one-loop (and only) self-energy diagram in XEFT. The diagram with external legs amputated can be expressed 
as -iS(F;cn,)(5*T 


where S is the mass difference in Eq. (2) and is the Dtt reduced mass defined in Eq. (3b). The branching fraction 
for D*^ —>■ D^tt^ has been measured to be (62 ± 3)%, but the full width of has not been measured. 

By exploiting isospin symmetry, the coupling constant can also be determined from decays of the D*~^. 

The full width of has been recently measured with high precision by the Babar Collaboration [17]: E[£)*+] = 
83.3 ± 1.2 ± 1.4 keV. The sum of the branching factions of into Zl°7r+ and Zl+7r° is 0.984 ± 0.005. In Galilean- 
invariant XEET, the sum of the partial widths into Zl°7r+ and is 


r[D*+ D°Tr+] + T[D*+ D+Tr°] 


1 

37r 



2(2fx^S+of/^ + (2/x^(5++)^/^ , 


(19) 


where (5+o = 5.856 ± 0.002 MeV is the D*^ — I?°7r+ mass difference, and <5++ = 5.68 ± 0.08 MeV is the D*'^ — 
mass difference. The resulting value of the coupling constant g is given by 

a‘2 

= (3.67 ± 0.08) X 10"® MeV^ (20) 

4m/2 


The pion decay constant determines the scattering amplitude for tttt scattering. With the convention for used 
in Ref. [6], its value is /tt ~ 132 MeV. Given that value, the pion transition coupling constant in Eq. (20) is determined 
to be « 0.59. A more appropriate dimensionless measure of the strength of the pion exchange interaction in XEET 
is 


127rTO/2 


0.15, 


( 21 ) 


where /i is the reduced kinetic mass of D*^D^ defined in Eq. (3a). 

The determination of the values of g and separately is unnecessary in XEFT applied to the DDtt sector (which 
includes D*D), because the number of pions in the system can only be 0 or 1. In the DDtttt sector (which includes 
D*D*), the number of pions in the system can be 0, 1, or 2. Since there can be contributions from tttt scattering, the 
value of /tt is needed at a sufficiently high order in the XEFT power counting. 


B. Propagator 

The conservation of pion number in Galilean-invariant XEFT implies that the exact D*^ propagator can be cal¬ 
culated analytically. It can be obtained by summing a geometric series of the one-loop D*'^ self-energy diagram in 
Figure 5. Galilean invariance implies that the propagator for a with energy pq and momentum p is a function of 
the Galilean-invariant combination 


Acm — PO 

The exact D*^ propagator can be expressed as 

Dl^{po,p) = 


2{M + m) 

i[l + 6Z]6^^ 


Ecm — (E* -|- SE^:) — S(ii/cm) 


( 22 ) 


(23) 


where S(Ecm) is the self-energy and is a rest-energy counterterm. We have introduced a wavefunction 
renormalization factor 1 -|- 6Z in the numerator. It can be made completely arbitrary by absorbing canceling factors 
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~ >- 


FIG. 6: The self-energy counterterm vertex. The Feynman rule for this vertex at order in the complex on-shell scheme 
is given in Eq. (31). 


of [1 -I- SZ] into the interaction vertices. With dimensional regularization, the D*^ self-energy is 


SiE, 


— Mtt ‘ 


\^mfl 


A3-^r(-d/2) „ 

(4^)d/2 [e 


,d/2 


(24) 


where A is the renormalization scale and is the reduced mass in Eq. (3b). The self-energy has cubic and 

linear nltraviolet divergences that are manifested as single poles in d and d — 2, respectively. The pole in d — 2 is 


1 

d-2 


lim(d - 2)T,{Ecm) 


9 ^ \ 2/r^A p 

4m/2 7 37r(d-2) 


(25) 


In the physical dimension d = 3, the self-energy is pure imaginary for real positive Ecm- 


lim E,{Ecm) 



Mtt 

Ott 




(26) 


Its value at E = d is —fr*o,7r/2, where F^o.tt is the partial width of D*° into given in Eq. (18). 

In the minimal power-divergence subtraction (PDS) renormalization scheme [16], the rest energy and the 
propagator counterterms are 


A'*,PDS 
dZpDS 


<5, 

0 , 

( 9^ \ 2m^A 
\4mf^J 37r(d- 2)' 


(27a) 

(27b) 

(27c) 


In the physical dimension d = 3, the self-energy S(Ecni) reduces to Eq. (26). The approximate position of the pole 
in po of the D*^ propagator in Eq. (23) is where E^m is equal to the complex energy 5 — fr*o,7r/2. The PDS scheme 
does not take into account the decay —>• 

The pole in the energy po for the physical propagator occurs when E^m is equal to the complex energy 
E^: = S — ir*o/2, where r*o is the full width of the id*°. We introduce the complex on-shell (COS) renormalization 
scheme for the D*'^ propagator in which its pole in po is at the complex physical value and the residue of that pole is 
the same as at LO. The rest energy and the propagator counterterms are 


E^: = d —fP^o/S, (28a) 

SE^ = -E{E^), (28b) 

dZ = -E'(E*), (28c) 


where E(i?) is the self-energy in Eq. (24). In the physical dimension d = 3, dE* and 6Z in the COS scheme reduce to 

If we ignore the tiny difference between E^, = 6 — ir*o/2 and d, the rest energy counterterm in Eq. (29a) reduces to 
8E^ = —ir*o,-n-/2, where P^o.tt is the partial width of the into in Eq. (18). 

Although Eq. (23) is the exact propagator in Galilean-invariant XEFT, it is not practical to use this propagator 
in loop diagrams, because the term in the denominator proportional to E^!^ makes loop integrals more complicated 
to evaluate. It is better to expand the exact propagator in Eq. (23) in powers of 


D"J{Po,p) 


i6^3 

Ecm 


oo 


E 

N=0 


Z^iEcm) + dE* -|- 6Z{Ecm ~ A"*) 
[l + SZ]{Ecm-E,) 


(30) 



11 



FIG. 7: The one-loop diagram for the LO D*^D^ transition amplitude. There is an implied sum over the two possible directions 
of the arrows on the internal lines. The diagram with external legs amputated can be expressed as the contact interaction 
vertex (—1(70/2)5*^ multiplied by (7011(170111), where n(17cm) is the LO D*^D^ self-energy. 


The expansion generates a geometric series of propagator corrections. The terms of order correspond to inserting 
into the propagator the one-loop D*^ self-energy subdiagram in Fig. 6 and a D*^ self-energy counterterm. The 
Feynman rule for the order-g^ self-energy counterterm vertex in the COS scheme is 


+ ^’{E,)[po-p‘^/{2{M+ m))- 


(31) 


where po and p are the energy and momentum of the D*^, i and j are its vector indices, and E(E) is the self-energy 
in Eq. (24). 

The previous NLO calculations in XEFT were carried out in the original version of XEFT, which is not Galilean 
invariant. The D*^ self-energy E(po,t) is therefore a function of the two independent variables po and p. The kinetic 
mass of the D*^ was set equal to its physical mass M*. The rest energy was taken to be the real energy <5, 
although it was actually taken into account through the rest energy term for the 7r°. In the NLO calculation of 
the decay of X(3872) into in Ref. [6], the renormalization of the propagator was carried out using the 

minimal power divergence subtraction scheme in Eq. (27). In the physical dimension c? = 3, the counterterm (5 Zpds 
in Eq. (27c) reduces to the renormalization scale A multiplied by a constant. This term was cancelled by a kinetic 
counterterm for the H*°. 

In the NLO calculation of the D*^D^ scattering length In Ref. [11], the complex energy of the 77(3872) and its partial 
width into were calculated to NLO using original XEET. The calculation involved a two-loop diagram that 

has a one-loop 77*° self-energy subdiagram and that diverges as 1/p as the relative momentum p of i7*°i7° approaches 
0. The divergence was avoided by resumming to all orders terms proportional to F^^q^^ from one-loop 77*° self-energy 
subdiagrams. The resummation could have been implemented in such a way that the terms that were resummed had 
only integer powers of F^q.tt- This is equivalent to the complex on-shell renormalization scheme for the 77*° propagator 
in Eq. (28), but with 77* replaced by i5 —zr*o,,r/2. The imaginary part of 77* is not equal to the physical value —F*o/2, 
because the decay 77*° —>■ 77°7 was not taken into account in Ref. [11]. 


C. D*°D° Transition Amplitude 

In XEFT, the conservation of pion number ensures that loop diagrams in the DDir sector have a single 77 or 77 
propagator in every loop. The integral over the loop energy of a 77 can therefore be evaluated by closing the contour 
around the pole in the 77 propagator, putting the 77 on its energy shell. Since 77 or 77 lines in loops can be put on 
their energy shells, it is sufficient to consider Green functions in which all external 77 or 77 lines are on their energy 
shells. 

The amputated connected Green function for 77*°77° —>■ 77*°li° in the C = + channel, with the incoming and 
outgoing 77° or 77° on their energy shells, can be expressed as , where the tensor is a function of the total 

energy Pq and the momenta of the two incoming and two outgoing particles. Its vector indices are associated with 
the polarizations of the incoming and outgoing 77*° or 77*°. In XEFT, the corresponding amplitude in the C = — 
channel is 0. We will refer to A^^ as the transition amplitude for D*'^D^ in the C = + channel or simply as the 
transition amplitude. Its contribution to the amplitude for 77*°77° —>• 77*°77° is -|-z-4*'’72, where the factor of (l/-\/2)^ 
comes from projections onto the C = + channel. 

The tree level term in the LO transition amplitude, which is given by the Feynman rule in Eq. (13), is —CqS^E 
The 77*°77° coupling constant Co must be treated nonperturbatively in order to generate the bound state that can 
be identified with the A(3872). The one-loop diagram in Fig. 7 is therefore also LO. The complete LO transition 
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FIG. 8: The Lippmann-Schwinger integral equation for the LO D*^D^ transition amplitude. In the second diagram on the 
right side, there is an implied sum over the two possible directions of the arrows on the internal lines. The Feynman rule for 
the LO transition amplitude is given in Eq. (44). 


amplitude can be obtained by summing a geometric series of the one-loop diagrams. This sum is equivalent to the 
solution of the Lippmann-Schwinger equation shown in Figure 8 . 

In Galilean-invariant XEFT, the LO transition amplitude depends on the energies and momenta of the initial-state 
and final-state particles only through the Galilean-invariant combination 


Erm — — 


pi 


2(2M + mY 


(32) 


where Pq is the total energy of the D*^ and relative to the threshold and P is their total momentum. 

The LO transition amplitude is (2.7:/ij)A{Ecm)Y where 


A{Ecm) — 




-[CY^ +SCY']+n{E^m) 


(33) 


and n(i?cm) is a function of Ecm that we will refer to as the D*^D^ self-energy. The inverse coupling constant from 
the contact interaction vertex in Eq. (13) has been separated into and a LO counterterm ■ Despite the 
notation, 60^^ is independent of the NLO counterterm SCq in the counterterm vertex in Eq. (16). With dimensional 
regularization, the D*^D^ self-energy is 


U{Ecm) = - Ecm)f"~\ 


(34) 


where is the reduced kinetic mass of D*^D^ defined in Eq. (3a). If the amplitude A{Ecm) in Eq. (33) has a pole in 
P 07 the residue of the pole is 


(4^)^/2-iArf-3 

^ 2 /xr (2 - d/ 2 ) 

This is a smooth function of d for d < 4. 

The D*^I)^ self-energy has a linear ultraviolet divergence that is manifested by a pole in d — 2: 

—^ lim (d - 2)n(Ecm) = 
d-2d^2 7r(d-2) 

In the minimal power-divergence subtraction (PDS) scheme [16], this pole is cancelled by a LO counterterm: 

The LO transition amplitude in Eq. (33) then has a finite limit as d —>■ 2: 

\imA{Ecm) = log (2yi{E,-Ecm)/ADA' 

The momentum scale A 2 in the logarithm is determined by the limiting behavior of as d —>■ 2. 

In the physical dimension d = 3, the LO transition amplitude in Eq. (33) reduces to the form 

lim A(Ecm) = - , , 

<^->■3 —Jx + a/— 2^ (Ecm — EY 


(35) 


(36) 


(37) 


(38) 


( 39 ) 
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where is a complex constant that we will refer to as the binding momentum of the X(3872). It is determined by 
the limiting behavior of as d —>■ 3: 

= (27r/M) hm + 6C^^) . (40) 

The amplitude A^Ecm) in Eq. (33) has a pole in Pq on the physical sheet. The residue of the pole is 

lim Zx =-7x//i- (41) 


The pole in Pq of the amplitude A(Ecm) in Eq. (33) occurs when Ecm equals a complex energy Ex that can be 
interpreted as the energy of the X(3872) resonance: 

Ex = E,-^\l{2ii) (42a) 

= (5-(5x-frx/2. (42b) 

In the first line, the complex energy Ex has been expressed in terms of the complex D*'^ energy E* and the complex 
binding momentum 7x. In the second line, Ex has been expressed in terms of real variables bx and Tx that can 
be interpreted as the binding energy and full width of the X(3872). The measured value of the binding energy bx is 
given in Eq. (1). The width Ex has not yet been measured, but there is an upper bound: Fx < 1-2 MeV [18]. If the 
complex binding momentum is expressed as 7x = 7re + *7imj where 7re and yim are real, the binding energy and full 
width are 

bx = (7re-7im)/(2Ai), (43a) 

Tx = r*o + 27re7im/Ai- (43b) 


In calculations beyond LO in the XEFT power counting, it is convenient to replace the contact interaction vertex 
in Eq. (13) by a vertex for the transition amplitude. The Feynman rules are the same for the D*^D^ D*^D^, 

D^D*^ O and D^D*^ —>• D^D*^ transition amplitudes: 


+i 



—AiE,^W^, 


(44) 


where E^m is the function of the total energy Pq and the total momentum P of the charm meson pair defined in 
Eq. (32) and i and j are the vector indices of the incoming and outgoing D*^. In d dimensions, the function A{Ecni) 
is defined by Eqs. (33) and (34). In the physical dimension d = 3, it reduces to Eq. (39). The vertex for the transition 
amplitude is represented by a blob, as illustrated for the D*^D^ —>■ channel by the Lippmann-Schwinger 

equation in Fig. 8. 

Beyond leading order in the XEFT power counting, the transition amplitude A^^{Ecm,P,p') for D*^D^ in the 
C = + channel in Galilean invariant XEFT is a function of Ecm and the relative momenta p and p' of the incoming 
and outgoing charm mesons. At NLO, this amplitude has additional ultraviolet divergences. The renormalizability 
of XEFT requires that they be cancelled by appropriate choices of the constants bC^ and bD^ in the counterterm 
vertex in Eq. (16). A renormalization prescription for the amplitude A^^ {Ecui,p,p') is necessary to determine the 
finite parts after the cancellations. In order to take into account decay modes of the A(3872) with momenta too large 
to be described explicitly, it is advantageous to use the complex on-shell (COS) renormalization scheme in which the 
pole in the total energy Pq is at its physical value. The complex energy of the pole can be specified in terms of the 
complex binding momentum yx of the A(3872), as in Eq. (42a), or in terms of its binding energy bx and full width 
Tx, as in Eq. (42b). The definition of the COS scheme for the D*^D^ —>■ D*^D^ Green function can be completed by 
specifying the residue of the pole in Pq. The amplitude A^^ {Ecni,P,p') has a well-behaved limit as p —>■ 0 and p' —>• 0 
that is a function of Ecm only. Thus we complete the definition of the COS scheme by requiring the residue of the 
pole in Pq for the amplitude (Acm, 0,0) to be the same as for the LO amplitude (27r/p)A(Ecm)^*A In the physical 
dimension d = 3, the limiting behavior as Pq approaches the pole is 


A^^{E,^,0,0) 


{2TTlp)Zxb^^ 


[Pq - PV(2(2M + m))] - [A, - 7i/(2p)] ’ 


(45) 


where Zx is the LO residue in Eq. (41) and yx is the complex binding momentum of the A(3872). The COS scheme 
can be implemented through specific choices of the counterterms bCg and bDo. 

In the previous NLO calculations in XEFT in Ref. [6] and [11], the NLO counterterm bCo was included, but bDQ 
was not. The counterterm bCQ could have been chosen so that the pole in Pq was at the physical point, as in Eq. (45). 
However without the counterterm bDQ, the prescription in Eq. (45) for the residue of the pole could not have been 
implemented. In fact, the residue would have been ultraviolet divergent. 
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IV. SCATTERING 

In this Section, we use Galilean-invariant XEFT to calculate the D*^D^ scattering length to NLO in the XEFT 
power counting. 


A. NLO Transition Amplitude 


In Galilean-invariant XEFT, the NLO transition amplitude for D*^D^ —>■ D*^D^ in the C = + channel in the 
center-of-momentum frame is a tensor {E,p,p') that depends on the total energy E and the relative momenta p 
and p' of the incoming and outgoing charm mesons. The LO transition amplitude is (2'k/ iJ,)A{E)S'^^, where A{E) in 
the physical dimension d = 3 is given in Eq. (39). The NLO diagrams for the transition amplitude are calculated 
in Appendix B. The 20 diagrams are labeled An, Bn, Cn, and Dn, where n is an integer. The NLO terms in the 
transition amplitude have a well-behaved limit as p —>■ 0 and p' —?► 0 that is diagonal in the vector indices i and j. 
The NLO transition amplitude can be expressed as 


O-TT / 

A^i{E,Q, 0) = —(A{E) + {g^p^/Amfl)[F{E)A\E) + G{E)A{E)\ 

+{C2/c^)h{e)A^{e) - 

/ 


The NLO correction terms either have a factor of A?{E) or a factor of A{E). The function E{E), which has dimensions 
of momentum, is obtained by adding the pion-exchange diagram AA in Eq. (B8) and the D* propagator correction 
diagrams B\ and B2 in Eqs. (B9) and (BIO). The function B{E), which is dimensionless, is the sum of the pion- 
exchange diagrams A2 and A3 in Eq. (B7). The function E[(E), which has dimensions of (momentum)^, comes from 
the vertex diagrams C in Eqs. (BIG). The SCq and SDq terms come from the D*D counterterm diagrams D in 
Eq. (BIS). The functions F{E), G{E), and F[{E) can be expressed in terms of the loop integrals In, and Kimn 
defined in Appendix A: 


F{E) 


G{E) 

H{E) 


(7^= - 2p(2E, - rE)Kni (E) - (2 - r) Ji (E)^] 

+2[KiwiE) - 2pE,Ki2o{E)] + r/i(E,) [4pE* J2(E) - dJi(E)]), 

SrVr^liE - E,)Ji{E) - rEh{E)] 

(i[(I - r)E - E*] 

8TTg{E — Etf)Ji{E), 


(47a) 

(47b) 

(47c) 


where r = /i^/p is the ratio of reduced masses defined in Eq. (4). The terms proportional to A{E) in Eq. (46) have a 
single pole in E at E^: — 7^/(2p). The terms proportional to A^(E) have an unphysical double pole in E. The N^LO 
term would have a triple pole and higher order terms would have even higher poles. These unphysical multiple poles 
can be summed to all orders, in which case they produce a shift in the position of the pole in the LO amplitude. In 
Ref. [II], such a resummation was used in the NLO calculations of the binding energy of the A(3872) and its partial 
width into . It was not used in the NLO calculation of the D*^D^ scattering length. An expression for the 

transition amplitude that is accurate to NLO and has only a single pole is 


A*^(E,0,0) 


(27r/p) [I -h {g‘^p^/Amfl)G{E)] 

A(E)-i - {g^p^lAmn)F(E) - (C2lC^)iI(E) + (2i:jpGl)\8C^ + 8D^E \' 


(48) 


To NLO accuracy, the factor of I -I- {g‘^g?/Amf^)G in the numerator could equally well be moved to the denominator 
as a factor I — {g‘^p?/Amf‘^)G multiplying A{E)~^. 


B. Renormalization 

Loop integrals in XEFT have ultraviolet (UV) divergences. With dimensional regularization, the UV divergences 
produce poles in d — 2 and poles in d — 3. A pole in d — 3 represents a logarithmic UV divergence, and a pole in 
d — 2 represents a linear UV divergence. In previous NLO calculations in XEFT [6, II], power divergence subtraction 
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was used to remove the poles in d — 2. The subsequent limit d —>■ 3 produces terms that depend explicitly on the 
renormalization scale A. In the PDS scheme, a loop integral with a pole term (d — 2) requires a canceling 

counterterm —A/{d— 2). In the physical dimension d = 3, the loop integral becomes independent of A, but the 
canceling counterterm reduces to —A. 

In the NLO calculation of the decay of the X(3872) into in Ref. [6], the explicit dependence on A was 

through terms proportional to A, A^, and log A. The explicit dependence on A was exactly cancelled by the implicit 
dependence on A of the coupling constants of XEFT. In the NLO calculations of Ref. [11], in addition to terms 
proportional to A, A^, log A, there was explicit dependence on A from terms proportional to 1/(A — 7), where 7 is 
the leading order binding momentum or inverse scattering length. These terms were introduced by a resummation 
prescription for dealing with an infrared divergence at the D*^D^ threshold. The terms proportional to 1/(A — 7) 
were not cancelled by the implicit dependence on A of the coupling constants, but they were suppressed by a factor 
proportional to the partial width for —>■ I?°7r°. 

The poles in d — 3 in the NLO transition amplitude are given in Appendix B in Eqs. (B22). In the NLO transition 
amplitude in Eq. (46), the only poles in d — 3 come from the two-loop integrals ATno, K 120 , and Km in the function 
F{E) in Eq. (47a). The poles in Ki 2 [) and Km are constants given in Eqs. (A20b) and (A20c), respectively. The 
pole in ATiio is a linear function of E given in Eq. (A23). Thus the dependence on the energy in the pole terms is of 
the form J\?{E) and EJ\?{E). These poles in d — 3 can be cancelled by the dCo and 5Dq counterterms in Eq. (46). 
Thus logarithmic UV divergences in the NLO transition amplitude can be cancelled by the available counterterms. 

The poles in d — 2 in the NLO transition amplitude are given in Appendix B in Eqs. (B20) and (B21). In the NLO 
transition amplitude in Eq. (46), the functions G{E) and H{E) have single poles in d — 2 while the function F{E) has 
double poles. The single poles in d — 2 in the one-loop integrals Ji and I\ are given in Eqs. (A14). The poles in d — 2 
in the two-loop integrals are given in Eqs. (A18): ATno has a double pole and 1^120 has a single pole. The pole terms 
with the energy dependence J\?{E) and EA^{E) can be cancelled by the counterterms dCo and respectively. 
They include the double-pole terms from F{E), some of the single-pole terms from F{E), and the single-pole terms 
from H{E). However there are also single-pole terms from F{E) that have the energy dependence log(E* — E)A{E)‘^. 
In the single pole of G{E) in Eq. (47b), the energy-dependent denominator is cancelled, so it gives pole terms with 
the energy dependence A{E). Neither of these terms can be cancelled by the available counterterms SGq and 5D^). 
The poles in d — 2 after suitable choices of the counterterms are 


A^^{E) 


/ ^ 4rVl - r^A^ 2/x(E^. - E) 

V4m/2 ■ ' 

( 


\4.mfl 


d-2 
(—4)r-\/l — r g.A 
d-2 


A2 


-A^{E) 


A{E) 


(49) 


The argument of the logarithm has been made dimensionless by using the renormalization scale A. 

The uncanceled poles in d — 2 in Eq. (49) seem to indicate that the NLO transition amplitude in Galilean-invariant 
XEFT has linear UV divergences that cannot be removed by renormalization. This puzzle can be resolved by taking 
into account the expression for the amplitude A{E) in d = 2, which is given in Eq. (38). Its reciprocal A{E)~^ is 
linear in log(E* — E). The cancellation between the two terms in Eq. (49) leaves a term with the energy dependence 
A^{E) that can be cancelled by the counterterm 5Co- 

We proceed to implement the complex on-shell renormalization scheme specified by Eq. (45). This scheme requires 
the pole in E of A^^{E, 0, 0) to be at the same complex energy Ex and to have the same residue as the LO amplitude 
{271 /g)A{E)5^F It can be implemented as specific choices of the counterterms 5Gq and 5Dq. We use the variation 
of the amplitude in Eq. (48) in which the term {g^g?/Amf^)G in the numerator is moved to the denominator. The 
renormalized expression for the transition amplitude is 


^*^(E,0,0) = {27T/gi)F^[A{E)-^ - {g^g?/Amfl)[G{E)A{E)-^ - G{Ex){E - Ex)/Zx] 

-{g^gyAmfl)E,^^{E) - (C2/Co)idsub(i?)] (50) 

where Fsuh{E) and idsub(d?) are obtained by subtracting terms from F{E) and H{E): 

FsnbiE) = F{E) - E{Ex) - F'{Ex){E - Ex), (51a) 

id,ub(£l) = H{E) - H{Ex) - H'{Ex){E - Ex). (51b) 


The denominator in Eq. (50) vanishes ai E = Ex- The value of its first derivative at E = Ex is such that the residue 
of the pole in E is {277/g)Zx5‘^G in accord with the renormalization prescription in Eq. (45). The explicit expressions 
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for the counterterms are given by 


27rdGo 

_ pV" 

pGl 

1 

27rdL»o 

_ pv" 

pGl 



(^F{Ex) — \F'{Ex) + G{Ex)/Zx'\Ex'^ 
E^^(h{Ex)-H'{Ex)Ex), 
[E'{Ex) + G{Ex)/Zx] + ^H'{Ex). 


(52a) 

(52b) 


In the previous NLO calculations in XEFT in Refs. [6] and [11], power divergence subtraction was used to make 
double and single poles in d — 2 explicit as terms proportional to and A, where A is the renormalization scale. 
Since the poles in d — 2 have the same energy dependence as the counterterms 5Co and dHo, the terms proportional 
to A^ and A are exactly cancelled by those counterterms in the COS scheme. Thus it is unnecessary to use power 
divergence subtraction to make the poles in d — 2 explicit. 


C. Scattering Length 

T-matrix elements in XEFT involving charm mesons and 7r*^’s can be obtained from the amputated connected Green 
functions by the following steps: 

• Multiply by a polarization vector e* for a or in the initial state and by a polarization vector for a 
D*^ in the Hnal state. 

• Multiply by a residue factor for each D*^ or D*^ in the initial and final states. In the COS scheme, 

Z, = 1. 

• Put the external lines on their energy shells. For the or and for the 7r°, the energy shells are real: 
Po = P^/(2M) and po = P^/(2 to), respectively. For the or Z)*°, the energy shell is complex: po = 
E^, +p^/(2(M + to)). 

We consider the elastic scattering of D*^D^ in the center-of-momentum frame with incoming momenta Epz and 
outgoing momenta ±pn. The scattering angle satisfies cos 9 = z ■ fi. The initial and final polarization vectors of the 
D*^ are e and e'. The total energy is Ep = F* +p^/(2p). The T-matrix element is 

T{p,e) = {l/2)Y,A^^{Ep,pz,pn)Fe'^*. (53) 

'ij 

In the limit of zero relative momentum, the scattering is isotropic. The T-matrix element reduces to 

Tip = 0) = (I/2)^M*^(F;*,0,0)eV^*, (54) 

ij 

where the transition amplitude is given by Eq. (50) evaluated at E = E^,. 

In the case of short-range interactions, a scattering amplitude can be expanded in powers of the relative momentum. 
This expansion is called the effective range expansion. Beyond LO in XEFT, the effective range expansion for D*^D^ 
breaks down because of the effects of the exchange of a pion that can be on its energy shell. However Jansen, Hammer, 
and Jia pointed out that the leading term in the effective range expansion, which is the S-wave scattering 

length, remains well defined [11]. They calculated the scattering length to NLO in original XEFT, truncating the 
expression at first order in an expansion in powers of yjc / V 2m5 and at leading order in m/M. We will calculate the 
scattering length to NLO in Galilean-invariant XEFT, truncating the expression at fourth order in an expansion in 
powers of 7x/-\/2pi?* but keeping all orders in m/M. 

The S-wave scattering length 0 ^,+ for D*^D^ in the G = + channel can be defined by expressing the T-matrix in 
the zero-momentum limit in Eq. (54) as 

r{p = 0) = (l/2)(27r/p)(-a«.+)£ • e'*. (55) 

At leading order in XEFT, the inverse scattering length is equal to the complex binding momentum: l/as^+ = 7x- 
The inverse scattering length at NLO is 

1K,+ = 7x + (5V/4m/2)[G(E;x)-2G(E;,)]7x/2 

+ (pVV4m/2)F,,b(i^*) + (G2/Go)i/.ub(£;*). 


(56) 
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We have used the expressions for Ex and Zx in terms of jx in Eqs. (42a) and (41). This expression for \/as,+ in 
Eq. (56) can be simplified by expanding it in powers of qx/VS/xE* using the expressions for the one-loop integrals J„ 
and In in Eqs. (A21a) and (A22a) and the threshold expansions for the two-loop integrals Kimn in Eqs. (A30). The 
expansion of G{Ex) — 2G{E^,) to third order in yx is 

GiE,) - 2G(E.) « -|i,.=/=^rV7y2^(l + , (57) 

where r = is the ratio of reduced masses defined in Eq. (4). The quantity Hsuh{E*) is very simple: 

i7sub(A,) =7^/2. (58) 


The expansion of Esub(E'*) to fourth order in yx is 

J.5/2J1 _ _ j.] 


Esub(A*) 


-l- 


QTTy/l — r 


7x a/ 2/xA* 


12x2 


j,5/2 ^ (12 — 22r — 5r2)r^/2 ^ (2 — r)(2 — 4r — r2) arccos(-yr) 


8(2-r) 
\/l — r 




7x 

2^E^ 


(59) 


Our final result for the inverse scattering length, including all terms suppressed by up to three powers of yx/•\/2/ii?*, 
is 


l/as,+ = 


G2 




lx + :^7=r7x ' n 

2Go 12xm/2 


A/2 fi_ .^J,2^ 7x 


8 J 


--r- -r^ _Lr5/2 _lx 


1- r - r~ - r' 

2 8 lOx 


4 

2^Et 


(60) 


The coefficient of each term has been expanded to relative order . The leading power of r in the g^yx term comes 
from the function G. The entire g^lx term comes from the function G. The <7^yx term receives contributions of 
order from both the function E and the function G. 


D. Comparison with Previous Calculation 

We can compare our result for the scattering length 0^^+ in Eq. (60) with that from a previous calculation by Jansen, 
Hammer, and Jia [11]. They used original XEFT to calculate the complex energy Ex of the A(3872) and the complex 
scattering length as,-i-: dropping terms that were suppressed by at least one power of r = m/M. They expressed their 
results in terms of a leading-order inverse scattering length y. Their physical renormalization condition was that 
the real part of the complex energy Ex should be identified with the energy of the A(3872) resonance. Their NLO 
corrections included terms proportional to g^, G 2 , and an additional coupling constant D 2 that takes into account 
the dependence of the contact interaction on the pion mass. The coupling constant D 2 allowed them to study the 
dependence of the mass of the X (3872) on the light quark masses. Their NLO corrections also included terms involving 
the partial width for —>■ that were resummed to all orders. If the other NLO corrections are omitted, 

their results reduce to 

Ex = J-yV(2A)-/r*o.^/2, (61a) 

as,+ = 1/ (7 - , (61b) 

where r»o,ir is the partial width of D*^ into The universal relation expressing the binding energy in terms of 

a large positive scattering length is 


Ex = S-l/{2fial^). 


( 62 ) 
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The leading order results in Eqs. (61) are incompatible with the analytic continuation of this universal relation to a 
complex scattering length as^+. This incompatibility suggests that the resummation used in Ref. [ 11 ] to eliminate an 
infrared divergence at the D*^D^ threshold was inadequate. Because of this problem, we will only compare results in 
the limit ^ 0. In this limit, their NLO results reduce to 


Ex = 


as,+ = 




1 + 


(C'2/C'^)(47r7/^) + (52^/67r/2)(l/7)[(A - - 2m(I4] 


1 

7 


2/x 

1 f 


Vl27r/2 


1 + {g'^fi/6Trfl^){mS/j) log(l + 2i^jx'2m5) 
—2i^\/2mS + — 2jA — 2m5d'2 


(63a) 

(63b) 


where A is the PDS renormalization scale and d '2 is a finite constant related to the coupling constant 02- Note that 
the dependence on A can be eliminated from both Ex and from as,+ by absorbing it into In. the expression for 
as,+ in Eq. (63b), the factor of 1 / 7 ^ comes from a term that has an unphysical double pole in the energy. The double 
pole could have been eliminated in favor of a shift in the position of the single pole by expressing as^+ as the ratio of 
a numerator and denominator with NLO accuracy, as in their result for Ex in Eq. (63a). 

The complex on-shell renormalization scheme can be implemented by choosing d '2 so that the numerator of the 
NLO correction term to Ex in Eq. (63a) is zero. The complex energy Ex then reduces to d — 7 ^/( 2 ^), so 7 can be 
identified with the binding momentum 7 x of the AT(3872). The NLO scattering length calculated in Ref. [11] reduces 
to 


^s,+ — 

lx 


2irC2 

+ 77 


9 d 
127r/^ 


(^2iV2md 


lx 


(64) 


This can be compared with the result from Galilean-invariant XEFT in Eq. (60). Dropping terms suppressed by three 
or more powers of jx / ^/2jlE^ and keeping only the leading power of r in the coefficients, our result reduces to 


as.+ = — 
lx 


1 - 


^ 2 , 

2Co^^ 127rm/2 




lx 


^/2gIE^ 


(65) 


If we ignore the tiny difference between E* and 5 and if we approximate /i by Af/2 and by m, the •\/2/ii?* term 
in Eq. (65) agrees with the ^/2m5 term in Eq. (64). However the other correction terms in Eq. (64) are proportional 
to 7 x, while those in Eq. (65) are proportional to Furthermore the coupling constants C 2 and Cq appear in the 
combination C^jC^ in Eq. (64), while they appear in the combination C^jC^ in Eq. (65). Thus there seems to be a 
clear discrepancy between the results for 0 ^^+ in Eqs. (64) and (65). However the conclusion that there is a discrepancy 
relies on an implicit assumption that the coupling constant C 2 in Eqs. (64) and (65) can be identified. In the following 
section, we will show the surprising result that in the limit g —>■ 0 in which the pions decouple, the NLO correction 
terms proportional to C 2 in Eqs. (64) and (65) are indeed compatible. Whether the apparent discrepancy between 
the NLO correction terms proportional to g^ix in Eq. (64) and to g^l\ in Eq. (65) is real cannot be determined 
definitively using the results given in Ref. [11]. 


E. Decoupled Pions 


If pions are decoupled by setting <7 = 0, the NLO terms in the transition amplitude {E,p,p') reduces to the 
vertex diagrams C in Eqs. (B16) and the D*I) counterterm diagrams D in Eq. (B18). The terms proportional to 
A^{E) which have a double pole in E, can to NLO accuracy be replaced by terms in the denominator, as in Eq. (48). 
After implementing the COS renormalization scheme, the transition amplitude reduces to 


A'^{E,p,p') 


(27r^)[l -7 ^{C2/Co){p^ +p1] Hi 
A{E)-^ - {C2/Co)H,ub{E) 


( 66 ) 


where HsuhiE) is defined by Eqs. (51b) and (47c). 

To obtain the T-matrix element in Eq. (53), we put the charm mesons on their energy shells by setting p' = p and 
setting the energy equal to Ep = E* +p^/{2p). The term Hsuh{Ep) in the denominator of Eq. ( 66 ) is 

Hsnh{Ep) = \ix + llxP^ + ip^- (67) 


The T-matrix element reduces to 


r{p) 


_ (7r//i)[l -t (C 2 /Co)p^] _ 

-lx - ip - {C 2 /Co)[\'y\ + lixp"^ + ip^^ 


e 




( 68 ) 
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To NLO accuracy, the factor 1 + in the numerator of Eq. (68) could equally be included as a factor 

1 — (C' 2 /C'o)p^ multiplying —7x — in the denominator: 

-[7x + 5(C'2/C'o)7i] - 5(C'2 /C'o)7xP^ - ^ 

We can read off the S-wave inverse scattering length and effective range in the C = + channel as coefhcients in the 
denominator: 


l/ns,+ — 7x + 5(C'2 /Co)7x) (70a) 

r,,+ = -(C2 /Co)7x. (70b) 

The inverse scattering length in Eq. (70a) is consistent to NLO accuracy with setting g = 0 in Eq. (60). Eliminating 
C' 2 /C'o, we reproduce a familiar relation between the scattering length, effective range, and binding momentum to 
first order in the range expansion: 

l/as.+ = 7x (1 - ■ (71) 

Our result for the effective range in Eq. (70b) implies that, in the absence of pions, the ratio C' 2 /Co of coupling 
constants has a simple physical interpretation in terms of the effective range and the binding momentum: 


C 2 /C 0 = -rs,+ hx- (72) 

In the previous NLO calculations in XEET in Refs. [6] and [11], power divergence subtraction was used and Co and 
C 2 were determined as functions of the renormalization scale A. The dependence on A canceled in the following 
combination of coupling constants: 


C2/C^ = Mr«.+/(47r). (73) 

The apparent incompatible with our result in Eq. (72) is at first puzzling. Our coupling constants Cq and C 2 can 
be defined as the coefficients of the terms in the Lagrangian given in Eq. (9). This seems to be equivalent to the 
definitions of Co and C 2 in Ref. [6] up to corrections suppressed by powers of r = The resolution is that the 

absence of the counterterm 5Do in the Lagrangian for original XEET in Ref. [6] implies that the definitions are not 
equivalent. To determine the relation between the coupling constants, it is necessary to compare physical observables, 
such as the effective range By equating in Eqs. (72) and (73), we can infer that the combination C 2 /C 0 in 
Ref. [6] should actually be identified with our ratio C 2 /C 0 multiplied by —^7x/(47r). This identification brings the 
NLO correction terms to 0 ^^+ proportional to (C 2 /Co) 7 x in Eq. (64) and (C 2 /Co) 7 x in Eq. (65) into agreement. 


V. SUMMARY 

We have introduced a Galilean-invariant formulation of XEET to exploit the fact that mass is very nearly conserved 
in the transition D* —>■ Dx. Galilean invariance requires the kinetic mass of the D* to be equal to the sum of the 
kinetic masses M and m of the D and tt. One advantage of Galilean invariance is that an amplitude is the same in 
every Galilean frame. A more important advantage is that it strongly constrains the form of ultraviolet divergences. 
The operators in the NLO Lagrangian for Galilean-invariant XEET can be obtained by making minor modifications 
in the terms of the NLO Lagrangian for original XEET. NLO calculations in Galilean-invariant XEET are ultraviolet 
finite for arbitrary values of the masses M and m. In contrast, NLO calculations in original XEET are ultraviolet 
finite only if they are expanded in powers of m/M and truncated at a sufficiently low order. 

We also introduced the complex on-shell (COS) renormalization prescription that takes into account the effects of 
transitions to states that cannot be described explicitly in XEET. These transitions include the decay —>■ , 

which accounts for about 38% of the total width of the Il*°. They also include the decays of Ar(3872) to all final 
states other than . These other decay modes may account for a substantial fraction of the total width of the 

AT(3872). The COS prescription requires as input the total width of the and the total width of the AT(3872), 
which has not yet been measured. 

We illustrated the use of Galilean-invariant XEET by calculating the transition amplitude for D*^D^ —>■ D*^D^ 
to NLO in the XEET power counting. We used dimensional regularization in d spatial dimensions to regularize the 
ultraviolet divergences. The NLO Lagrangian includes counterterms that can cancel terms whose energy dependence 
has the form and EA'^{E), where A{E) is the LO transition amplitude. It was straightforward to verify that 

the logarithmic ultraviolet divergences, which appeared as poles in d — 3, could be cancelled by the counterterms. 
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The cancellation of linear ultraviolet divergences, which appeared as poles in d — 2, was not as straightforward. In 
addition to the poles that could be cancelled by the counterterms, there were additional poles in d — 2 whose energy 
dependence had the form A^iE) log(£l* — E) and A{E). They cancelled each other only upon using the specific form 
of the LO transition amplitude A{E) for d = 2. 

We used our NLO result for the transition amplitude to calculate the D*^D^ scattering length as_+ to NLO. Our 
result differs in several respects from a previous result for as,+ calculated using original XEFT in Ref. [11]. One 
important difference is in the effect of the D*^ width. Our COS renormalization prescription ensures that the effect 
of the D*^ width is in accord with the universal relation between the binding energy of an S-wave bound state near 
threshold and the scattering length of its constituents. The effect of the partial width of the into in 

Ref. [11] is incompatible with the universal relation. This suggests that the resummation used in Ref. [11] to deal 
with an infrared divergence at the D*^D^ threshold was inadequate. There is an apparent disagreement between our 
result and that of Ref. [11] in the NLO corrections to as^+ proportional to C 2 , but we showed that the difference is 
actually due to different definitions of C 2 - Our NLO correction to 0 ^,+ proportional to was expanded in powers of 
'jx where 'jx is the binding momentum of the X(3872). The leading term is proportional to g^V^/rnS, and 

it agrees with that in Ref. [11] in the limit m/M —>■ 0. The next term in our expansion is second order in 7 ^, while 
there is a term in Ref. [11] that is first order in yx. The origin of this discrepancy has not been determined. 

Our Galilean-invariant formulation of XEFT removes obstacles to precise calculations of the properties of the 
X(3872) resonance. Our COS renormalization prescription takes into account the 38% branching fraction of D*^ into 
and the significant branching fraction of X'(3872) into decay modes other than In original XEFT, 

there was a serious limitation on the precision from the need to truncate the expansion in m/M in order to avoid 
ultraviolet divergences. Our Galilean-invariant formulation makes this expansion unnecessary. The expansion may 
still be useful to obtain simple results, but it can be carried out to whatever order is required for the desired precision. 

XEFT was invented to allow the effects of pion transitions on the X(3872) resonance to be taken into account 
systematically. XEFT can also be used as a framework for quantifying the effects of physics beyond the DDtt sector. 
This physics include charmonium states such as the Xcii^P), which couples to D*D, and ip{3770), which couples to 
DD. It also includes tetraquark cc mesons and hybrid cc mesons. The improved precision of Galilean-invariant XEFT 
increases the motivation for using this effective field theory to quantify the effects of physics beyond the DDtt sector. 
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Appendix A: Loop Integrals 

In this Appendix, we evaluate the loop integrals that are required to calculate the NLO transition amplitude using 
dimensional regularization. We also determine their poles in d — 2 and d — 3, where d is the number of spatial 
dimensions. 


1. Reduction to momentum integrals 

The loop integrals have D*, D, and tt propagators. Every loop has a D propagator, and it is convenient to take 
its energy and momentum to be the loop integration variables. A D propagator has a denominator of the form 
po — P^/(2M) -I- ie. The integral over po can be evaluated by closing the contour around the pole at po = p^/(2M), 
which puts the D on its energy shell. With dimensional regularization in d spatial dimensions, the measure for a 
momentum integral is 


= X.3-d 


{2ttY 


(Al) 


where A is the renormalization scale. Tensor reduction can be used to reduce the loop momentum integrals to scalar 
integrals. Whenever possible, it is best to express the momentum-dependent terms in the numerator as a linear 
combination of the propagator denominators. After canceling terms in the numerator with propagators, the loop 
integrand reduces to a product of propagators. 
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The propagators in the loop integrand depend on the complex energy and on the masses m and M of the 
7r° and The loop integrals depend on the masses through the reduced D*D kinetic mass /i defined in Eq. (3a) 
and the reduced Dtt mass defined in Eq. (3b). It is convenient to express the dependence on the masses in terms 
of fj, and the ratio r = defined in Eq. (4). 

In loop integrals that have cuts with DDtt intermediate states, after integrating over the D and D energies, the 
pion propagator depends on the momenta p and q of the D and D and on the masses m and M. In evaluating the 
loop integral, it is often convenient to express the pion propagator in a form that depends on reduced masses instead 
of masses: 


E - {p + q)‘^/{2m) - {p^ + q^)/{2M) + ie 

i 


E -(q + VI - rpY/{2p^) - p^/{2p) + it 


(A2) 


Unlike the first expression for the pion propagator in Eq. (A2), the second expression is not manifestly invariant under 
the interchange p ^ q. 


2. One-loop momentum integrals 


The one-loop momentum integrals whose integrands are a D* propagator raised to an integer power have the form 

I 


ME) = j 

Jp 


[p-^ -2piE- EM' 


The analytic result for this integral is 


ME) = 


- d/2) 

(47r)<i/2 


[2p{E, - U)]^/^-". 


(A3) 


(A4) 


The branch cut in E is determined by the negative imaginary part of the complex energy E^. 

The one-loop momentum integrals whose integrands are a tt propagator raised to an integer power have the form 


ME) = [ 

Jn 


1 


- 2p^E - ieY 


The analytic result for this integral is 




(A5) 


(A6) 


where r = Pt^/ p is the ratio of the reduced masses defined in Eq. (4). The branch cut in E is specified by the it 
prescription in Eq. (A5). Since we choose the zero of energy to be at the threshold, we need only consider 

positive values for the real part of E. 

The one-loop momentum integrals whose integrands have a single tt propagator and a D* propagator raised to an 
integer power are 

ME,p) = ^ [^2 _ 2p{E - EM (P + g)V(2m) + M+q^)/i2M) - E - it 


The function Lq can be expressed in terms of the integral U given by Eq. (A6): 

Lo{E,p) = rh{E-pVM- 


(AS) 


This result can be obtained most easily by using the expression for the pion propagator on the right side of Eq. (A2), 
and then shifting the integration variable q. The function Li can be expressed as a Feynman parameter integral: 


EiiE,p) = r dx{x+{l-x)r)^ 

X [ — 2{x -f (I — x)r)p{E — (I — x)E^) + xp'^] 


(A9) 
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We will need the expansions of these functions to first order in : 


Lo{E,p) = rIi{E) 


1 - 


(d-2)p2 


ApE 


0{p% 


Li{E,p) = 


Ji{E)-hiE) 

2fj,A 


1 + 


r [4r(l — r)E + (4 — 4r — d) A]p^ 


2d/iA^ 


r{d - 2) [4(1 - r)E - dA] h {E)p‘^ 

Sdp^EA^ 


0{p% 


(AlOa) 


(AlOb) 


where A= (1 — r)A — A*. The expansion for Lq can be obtained most easily by expanding the expression in Eq. (A8) 
in powers of The expansion for Li can be obtained most easily by expanding the integrand in Eq. (A7) in powers 
of p and averaging over angles to reduce it to scalar integrals. The scalar integrals can be expressed in terms of the 
integrals J„ and /„ given by Eqs. (A4) and (A6). They can all be reduced algebraically to Ji and Ii. 


3. Two-loop momentum integrals 


The two-loop momentum integrals whose integrands have a tt propagator and one or two D* propagators raised to 
integer powers are 

KlmniE) = ' ' ^ 


[p^ - ME - E^n^ - 2p(E - A,)]" 

X__ (All) 

[(p + q)2/(2m) + (p2 + g2)/(2M)-A-ze]'' ^ ’ 

The momentum integrals with a tt propagator and a single D* propagator can be expressed as integrals over a single 
Feynman parameter: 


Kno(A) = M dx (1 - xMI^[2p{xE, - E^^ 

f dxx{l — x)~^^‘^[2fi{xE:^ — 

Jo 


KME) ^ /■ 


(47r)" 


(A12a) 

(A12b) 


where r = The momentum integrals with a tt propagator and two D* propagators can be expressed as an 

integral over two Feynman parameters: 

ifin(A) = f dww[2piwE, - E)M^ 

M) Jo 

X f dt [l — n; -I- rw^t{l — t)] . (A13) 

Jo 

This function depends on r through the prefactor r‘^/2 ^nd through the integral over t. For a given dimension d, 
Ariii(A) can be expanded in powers of r. The expansion is not straightforward, because there are negative powers of 
r that arise from the zc —>■ 1 endpoint region. 


4. Poles in d — 2 


In a dimensionally regularized loop integral, poles in d — 2 are associated with linear ultraviolet (UV) divergences 
in 3 spatial dimensions. An UV pole in d — 2 can arise from a factor of r(l — d/2) or r(2 — d) obtained by integrating 
over a loop momentum. In the two-loop momentum integrals A'no and Ari 2 o defined in Eqs. (A12), the integral over 
the Feynman parameter x also gives an UV pole in d — 2. 

The one-loop momentum integrals Ji and Ii given by Eqs. (A4) and (A6) have single poles in d — 2. We will need 
the pole and the constant term: 


Ji{E) - 

A 

\ 1 

1 

_1_ 

log- 

^ 27r 

[d-2 2 

log — 

h{E) - 
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1 
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_1_ 

(^log 
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2tt 

[d-2 2 


T 

2pE 


(A14a) 


(A14b) 
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The momentum scale A in the logarithms is 


A = ^/^A, 


(A15) 


where 7 is Euler’s constant. We will need the constant term in the integral J 2 given by Eq. (A4): 

A 


J2{E) 


87r/i(E* — E) 

We will also need the pole term in the integral Lq in Eq. (A8): 

rA 


LoiE,p) 


'27r(d- 2)' 


(A16) 


(A17) 


The two-loop momentum integral ATiio defined in Eq. (A12a) has a double and single pole in d — 2, while K 120 defined 
in Eq. (A12b) has only a single pole: 


KiME) 

Ki20{E) 
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2rA^ 


(47r)2 [(d-2)2 ' d-2 
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(47r)2(d-2)^(E;-E*)' 


—2p{E — E*) 1 
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(A18a) 

(A18b) 


5. Integrals at the D*°D° threshold 


If the one-loop integral Jn{E) given by Eq. (A4) is analytically continued to d = 3 and then evaluated at E = E^, 
it has an infrared divergence for n > 3/2. However if it is evaluated at E = E* and then analytically continued to 
d = 3, it vanishes because the integral has no momentum scale: 

Jn{E,) = 0, (A19) 


The two-loop integrals Kimn{E) defined in Eq. (All) can be evaluated analytically at the threshold E = E*. The 
integrals Ariio(E 4 ,), A'i 2 o(E*), and A'iii(E*) each has a single pole in d — 3. Their values near d = 3, including the 
pole in d — 3 and the finite term, are 


Kiio{E^) 


Ki2o{E^:) 


KniiE,) 
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(A20a) 

(A20b) 


(A20c) 


The momentum scale A in the logarithms is defined in Eq. (A15). The pole in Kiio{E^,) is the sum of an UV pole 
and an infrared pole in d — 3. 


6. Integrals near d = 3 


The one-loop integrals Ji and J 2 at d = 3 are 


= -^[2piE,-E)]^/^, 
J2{E) = -^[2p{E,-E)]-^/^. 


(A21a) 

(A21b) 
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The one-loop integrals Ii and I 2 d = 3 are 


I^{E) = ir ^/2 1 {A22a) 

Att 

E{E) = 1 [2/x£;]-i/ 2. (A22b) 


The two-loop momentum integrals Kimm{E) near d = 3 are expressed in terms of Feynman parameter integrals in 
Eqs. (A12) and (A13). In order to obtain expressions for these functions near d—3, it is necessary to make subtractions 
of the integrands to make them integrable at x = 1. The resulting expression for Kiiq(E), including the pole in d — 3 
and the finite term, is 


Kno{E) = 


47 . 3 / 2 ^ 

(47r)^ 


- 2 -b ^ logr -blog {E-2E^) 

1 

—2A, — - / dx (1 — — xA*) log 

2 Jo 


E — xA* 
E — Et, 


(A23) 


The pole with its residue evaluated at A = A* is also the UV pole in d — 3 of K{E^). The difference between the pole 
of K{E^) in Eq. (A20a) and this pole is the infrared pole in d — 3 of K{E^,). The UV poles in d — 3 of Ki 2 o{E) and 
Ariii(i?) are the same as the poles in their values at E^ given in Eqs. (A20b) and (A20c). The differences between 
their values at E in Eqs. (A12b) and (A12b) and their values at E^, are finite: 


Ki20{E) 


Ki2o{E^) 


4j.3/2 

(47r)3 


log 


E — E^ 

Etf 



x) log 


E — xE^ 
E — Etf 


K^^^{E) 


Kixi{E,) 


j.3/2 

(47r)3 



w 

•\/l — w(l — w -I- rw'^/A) 


log 


E — wE^ 

(1 — w)E^ 


(A24a) 

(A24b) 


7. Threshold expansions 

Threshold expansions for the two-loop integrals Kimn{E) near d = 3 cannot be obtained simply by expanding the 
expressions in Eqs. (A23), (A24a), and (A24b) in powers of A —A*, because this generates infrared divergences in the 
Feynman parameter integrals. The threshold expansion for A'iii(E) for general d can be obtained in a straightforward 
way by expanding the defining integral in Eq. (All) in powers of A — E*. The expansion through second order in 
E — E^ is 


Kni{E) = Kni{E,) + [K 2 iiiE,) + 2Ki2i{E,)]2fi{E-E,) 

+ [KsiiiE,) + 2K22 i{E,) + Ki 22{E,) + 2Ki3i{E,)] [2fi{E - E,)f 

+0{iE-E,f). (A25) 

The reason the threshold expansion of Kiii{E) is straightforward is that the only momentum region that contributes 
is where both loop momenta are of order 

In contrast, the integrals Ariio(i?) and Ki 2 o{E) have contributions from regions where one loop momentum is of 
order {^\E — These contributions have an expansion in non-integer powers of E — E^. By inserting the 

alternative expression for the pion propagator in Eq. (A2) into the defining integral for Kimo{E) in Eq. (All) and 
then shifting the momentum q, it can be expressed as 

Kirao{E) ^ ^ _ 2^{E - E*)]™ Jg - 2fj,-^Et + r{p^ - 2fi{E - A*)) - ief ' (A26) 

The expansion of the pion propagator in powers of — 2p{E — A,) produces a sum of products of one-loop integrals 
of the form Ii^k{E<,)Jm-k{E). With dimensional regularization, the terms with k >m are zero because the integral 
over p has no scale. The threshold expansion for Kimo{E) is the sum of the terms with the noninteger powers 
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{E — and the terms with integer powers that can be obtained by expanding the integrand of Eq. (A26) 

in powers of if — if*. The expansions of Kiiq{E) and Ki 2 o{E) through second order in if — if* are 

ifiio(E) = Kiw{E*) + rIi{E,)JiiE) + [K2ioiE,) + KME*)]2^iiE-E,) 

+ [KsioiE,) + K 22 o{E,) + KME,)] [2^(£; - E,)f + 0{{E - E,f), (A27a) 

Ki 2 o(E) = ifi 2 o(E*) + rh{E,)J 2 {E) - l 2 {E,).h{E) + [K 22 o{E,) + 2ifi3o(A*)] 2^x{E - A*) 

+ [^32o(£^*) + 2K23o{E,) + SKME*)] [2tJ.{E - A*)]2 + 0{{E - E,f). (A27b) 


The integer powers of if — if* in the threshold expansions in Eqs. (A25) and (A27) are determined by the two-loop 
integrals Kimn{E^:). By inserting the alternative expression for the pion propagator in Eq. (A2) into the defining 
integral for Kimn{E) in Eq. (All), its value at the threshold E = E^ reduces to 

Ip Iq [g^]”[(q -I- Vl - rp)2 -|- - 2r/rif* - ze]' ^ ^ 

The denominators of the integral over q can be combined by introducing a Feynman parameter. After evaluating the 
integral over q, the integral over p can also be evaluated analytically. The integral over the Feynman parameter gives 
a hypergeometric function; 


^ J Y ^ T{l + m + n-d)T{d/2-m)r{d/2-n) 

*’ ~ (47r)3 V4 ^) r(or2(d/2) 

X 2Fi{d/2 - m, d/2 - n, d/2, 1 - r) 2/xE*] . (A29) 

If n = 0, the hypergeometric function reduces to Inserting the expression for Kimn{E^d) in Eq.(A29) into 

Eqs. (A27a), (A27b), and (A25), the threshold expansions through second order in if — if* are 


iFiio(E) = iFno(E*)+rii(E*)Ji(E) 


2r3/2(2^E*) / I , 2^E* , 
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+ 82 Ei(-i,-i,f,I-r)) +0{{E-EY)- 

The function arccos('yr) in Eq. (A30c) has an expansion in odd powers of 


arccos 




lj.1/2 _ J_j,3/2 
TT 37r 


The hypergeometric function in Eq. (A30c) has an expansion in powers of r: 


37r 


2 ri[ 2 , r) ^ |i Y^24: 


I 1 


(A30a) 
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FIG. 9: Pion-exchange diagrams for —>■ D*^. They consist of the tree diagram yll, the one-loop diagrams A2 and 

A3 in which either the incoming or the outgoing charm mesons interact through a LO transition amplitude, and the two-loop 
diagram AA in which they both interact through LO transition amplitndes. The absence of arrows on the charm mesons lines 
in the loop for the diagram A4 implies a snm over the two possible directions of the arrows. 


Appendix B: Diagrams for NLO Transition Amplitude 


In this Appendix, we give results for the individual diagrams that contribute to the transition amplitude -|-fA*^/2 
for D*^D^ —>■ D*^D^ in the C = + channel, with the external legs amputated and the initial and final on their 
energy shells. In the center-of-momentum frame, the amplitude {E,p,p') is a function of the total energy E and 
the relative momenta p and p' of the incoming and outgoing charm mesons, respectively. The LO amplitude is the 
solution to the Lippmann-Schwinger integral equation shown in Figure 8. The amplitude at NLO can be expressed as 


A^^{E,p,p') = i27^/^i)AiE)6^^ 

n 


where the sum is over the NLO diagrams and the amplitude A{E) in the LO term is 


A{E) 


M/(27r) 

- 2fiJiiEy 


(Bl) 


(B2) 


The one-loop momentum integral Ji is given in Eq. (A4). In the physical dimension d = 3, the amplitude A(E) 
reduces to the expression in Eq. (39). The NLO diagrams for the transition amplitude can be organized into five sets 
of 4 diagrams labeled A, B, C, and D. 


1. Pion-exchange diagrams 


There are four diagrams that involve the emission of a pion by D*^ and its absorption by D^. The four pion- 
exchange diagrams, which are labelled Al, A2, A3, and A4, are shown in Eig. 9. The amplitude for the tree diagram 
A1 is a tensor in the vector indices i and j that depends on the relative momenta p' and p: 


■AJi{E,p,p') 



-2 


E — {p'^ + p' )/{2M) — {p + p'Y/{2m) + ie 


M 


M + m 


P + P 


M 


M + m 


P +P 


(B3) 


J 
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This can be expressed in a form that depends on the masses only though reduced masses: 


A^ME,P,P') 


\ (-2)(V1 -rp + p'YWl -rp' + pY 

4m/2y E-p‘^/{2^) - {p' + VI - rpy/{2p„) + ie’ 


(B4) 


where r = p. Unlike the expression in Eq. (B3), the expression in Eq. (B4) is not manifestly invariant under the 
simultaneous interchanges p ^ p' and i o j. The amplitude for the diagram A1 vanishes in the zero-momentum 
limit p,p' —>■ 0. 

The amplitudes for the one-loop diagrams A2 and AZ in Fig. 9 are tensors in the indices i and j that depend on 
the relative momenta p' and p, respectively. They can be decomposed into transverse and longitudinal components. 
The Feynman rule for the LO transition amplitude is given in Eq. (44). The amplitudes for the diagrams AZ and A2 
are 


a'Ue.p) = 






,A{E) 


\Amf^J {d- 1)V1 - rp'^' 

X { [ - (A(E,p) -b rp'^/p)Lo{E,p) - 2p{A{E,pf - 2rE^p‘^/p)Li{E,p) 

+r(A(E,p) +pyp)ME)] - pV7/) 

+{d- l)[A{E,p)Lo{E,p) +2pA{E,p){A{E,p) + rp"^/p)Li{E,p) 

-r{A{E,p) - (1 - r)p'^/p)Ji{E)\p^p>/p"^^, (B5a) 

A"i^{E,p') = A%{E,p ^ p'). (B5b) 

The one-loop integrals Ji, Lq, and Li are given by Eqs. (A4), (A8), and (A9), and A is a linear function of E and p^-. 

A{E,p) = {l-r)E-E,-p‘^/{2p). (B6) 


The expression for the amplitude A^^^{E,p) in Eq. (B5a) has terms proportional to s-iicl p^jP/p^ that 

are not analytic functions of the momentum vector p in the neighborhood of p = 0. The limit p —>■ 0 can be obtained 
by using the expansions for Lq{E,p) and Li{E,p) to first order in p^ given in Eqs. (AlO). The nonanalytic terms 
cancel in the limit p —>■ 0. The diagram A2 has the same zero-momentum limit. The sum of the amplitudes for the 
diagrams A2 and A3 in the zero-momentum limit is 


A 


ij 

A2+A3 


{E,0) 


9^ \ IGTrrVl - r p[{E - E^)Ji{E) - rEhjE)] 
4771/2^ d[{l-r)E-E,] > 


(B7) 


The one-loop integral Ii is given in Eq. (A6). 

The two-loop diagram A4 in Fig. 9 is the sum of a diagram in which a pion is emitted by D*^ and absorbed by 
and a diagram in which a pion is emitted by D*^ and absorbed by . The amplitude for this diagram is diagonal in 
the vector indices i and j and depends only on E-. 

= (^) [^K^io{E) - 2p(2E, - rE)K,^,{E) 

-i2-r)Ji{Ef]A^{E)5^A (B8) 

The two-loop integrals ATno and Km are given in Eqs. (A12a) and (A13). 


2. D* propagator correction diagrams 

There are two NLO diagrams involving a correction to a D* propagator. The two diagrams, which are labelled B1 
and B2, are shown in Figure 10. The diagram B1 has a D* self-energy subdiagram inserted into the D* propagator. 
It can be reduced to 


A%,{E) 



(—327r2)r/i 


[Kiio(E) - 2pE,Xi2o(E)] A^E) 


d 


(B9) 
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FIG. 10: D* propagator correction diagrams for D*^D^ —>■ D*°D^. The two-loop diagram B1 has a D* self-energy subdiagram 
inserted into the D* propagator. The one-loop diagram B2 has a D* self-energy counterterm inserted into the D* propagator. 

The two-loop integrals Khq and ^120 are given in Eqs. (A12). The diagram B2 has a D* propagator counterterm 
inserted into the D* propagator. The Feynman rule for the self-energy counterterm is given in Eq. (31). With the 
complex on-shell renormalization scheme for the D*^ propagator, the diagram B2 can be reduced to 

- dME)]A^E)S^E (BIO) 

The one-loop integrals Jn are given in Eq. (A4). 


3. vertex diagrams 


There are four NLO diagrams with a vertex. The four diagrams, which are labelled Cl, C2, (73, and (74, are 
shown in Fig. 11. The amplitude for the first diagram (71 is just the vertex in Eq. (15): 


A^i^{E,p,p') = (-C2/2) +/) 


(Bll) 


The factor of p^ + p'^ depends on the relative momenta of the incoming and outgoing charm mesons. If there is a 
contact interaction between the incoming charm mesons, p is replaced by a loop momentum k. The amplitude for the 
diagram (72 in which the incoming charm mesons interact through a LO transition amplitude can be obtained from 
Eq. (Bll) by multiplying it by (47r//i) Ji(E)7l(E) and replacing by the loop-integral-weighted average of 

A'S2iE,p) = -2TrC2Ji{E)A{E) (B12) 

The loop-integral-weighted average of is defined by 

s,pnE-E.-p/(2,.n 

' ' /»l/lE-£.-F/(2rt] ■ ' 

With dimensional regularization, the loop-integral-weighted average is very simple: 

(k^) =2fi{E- E,). (B14) 

The amplitude for the diagram (73 is obtained from Eq. (B12) by replacing p' by p. The amplitude for the diagram 
(74 is 


A^S4{E) = -16TT‘^C2Ji{EfA^{E){k‘^)d^E 


(B15) 


The sum of the four diagrams (71, (72, (73, and (74 has a multiplicative factor 1 -|- 47rJi7l that can be simplified by 
using the expression for Aq{E) in Eq. (B2): 


TtC2 


A^^{E,p,p') = ^ {p^+p'^)A{E) + IQME - EME)A\E) 

pL-o L 




(B16) 
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FIG. 11: The vertex diagrams for D*^D^ —>■ D*^D^. They consist of the vertex Cl, the one-loop diagrams C2 and 
(73 in which either the incoming or the outgoing charm mesons interact through a LO transition amplitude, and the two-loop 
diagram (74 in which they both interact through LO transition amplitudes. 


4. D*D counterterm diagrams 


There are four NLO diagrams that can be obtained from the diagrams in Figure 11 by replacing the vertex by 
a D*D counterterm. The four diagrams are labelled Dl, D2, D3, and DA. The amplitude for the diagram D1 is just 
the counterterm vertex in Eq. (16): 


A'j^^ = -[6Co+SDoE]d^C 


(B17) 


The amplitudes for the diagrams D2 and D3 each differs from this by a multiplicative factor 47r JiTl. The amplitude 
for the diagram DA differs by two such factors. The sum of the four diagrams Dl, D2, D3, and DA therefore has a 
factor [1 -|- 47rJi7l]^ that can be simplified by using the expression for Ao{E) in Eq. (B2): 




Att^ISCo + SDqE] 


A^{E)S^E 


(B18) 


The 5Co term in this expression can also be obtained from the LO transition amplitude in Eq. (B2) by replacing Cq ^ 
by — SCq/Cq and expanding the amplitude to first order in 5Cq. 


5. Complete NLO amplitude 

The complete NLO term in the transition amplitude {E,p,p') is the sum of (A) the pion-exchange diagrams 
in Eqs. (B3), (B5), and (B8), (B) the D* propagator insertion diagrams in Eqs. (B9) and (BIO), (C) the vertex 
diagrams in Eqs. (B16), and (D) the D*D counterterm diagrams in Eq. (B18). In the zero-momentum limit p,p' —f 0, 
the pion-exchange diagrams reduce to the sum of Eqs. (B7) and (B8) and the vertex diagrams in Eqs. (B16) reduce 
to the single term with the factor of A^iE). 


6. Poles in d — 2 

The momentum integrals with poles in d — 2 are given in Section A 4. In the NLO term in the transition amplitude 
A^^ {E,p,p'), all the terms with a double pole in d — 2 have the tensor structure . All the terms with a single pole 
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in d — 2 also have the tensor structure (5*^, with the exception of terms from the diagram ^3 with the tensor structure 
and terms from the diagram A2 with the tensor structure fp'^. Upon using the identities Lq = rJi and 

Li = 0 for the pole terms, the diagram Ai in Eq. (B5a) reduces to 

AUE,p) [5^^ +{d- 2)py/p2]. (B19) 


Thus the single pole in c? — 2 has the tensor structure (5*-' even at nonzero momentum. 

We now consider the terms with the tensor structure 5*^ in the zero-momentum limit. The pion-exchange diagrams 
A2 and Ai in Eq. (B7) and the vertex diagrams C in Eq. (B16) have single poles in d — 2: 


^g(E,0,0) 


/ \ (-4)rVl -r pK 


V4m/2 

87rC'2A 
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A{E) y, 


(B20a) 

(B20b) 


The two-loop pion-exchange diagram AA in Eq. (B8) and the sum of the D* propagator correction diagrams B1 and 
B2 in Eqs. (B9) and (BIO) have double and single poles in d — 2: 
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7. Poles in d — 3 


The only NLO diagrams for the transition amplitude that have poles in d — 3 are the two-loop pion-exchange 
diagram Ad in Eq. (B8) and the D* self-energy insertion diagram B1 in Eq. (B9). The momentum integrals with 
poles in d—3 are the two-loop integrals Kuo, K 120 , and Km. The poles are given in Eqs. (A23), (A20b), and (A20c), 
respectively. The pole terms are 
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